Driven particle in a cloud of mobile impurities 
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Abstract. The dynamics of a test particle interacting with diffusing impurities in 
one dimension is investigated analytically and numerically. In the absence of an applied 
external force, the dynamics of the particle can be characterized by a distribution of 
monotonic excursions Ax, which scales as a power law with an exponent ta x = 4/3. 
When the particle is driven at a slow constant velocity, there is again a power law 
distribution for the monotonic changes of the force Ai* 1 , which is characterized by a 
similar exponent ta_f = 4/3. These results can be understood from the theory of 
random walks. 
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1. Introduction 

The interaction of driven particles, flexible lines and membranes etc. with disorder is an 
important topic in condensed matter physics [H [21 [3] . Usually, this disorder is taken to 
be quenched, or frozen, such that its properties do not change within the relevant time 
scales. However, under certain conditions, this changes as in the case of the diffusion 
of solute atoms in metallic alloys [U [5] or oxygen vacancies in superconductors [6 J . The 
mobile impurities play an important role in the dynamics of such systems, as evidenced 
for example by the Portevin-Le Chatelier (PLC) effect in solid solutions [7J. There, 
within a certain range of temperatures and applied strain rates, the dynamic interaction 
of lattice dislocations and diffusing solute atoms result in phenomena such as negative 
strain rate sensitivity of the flow stress, giving rise to macroscopic serrations in the 
stress strain curve and strain localization in the form of bands of activity of various 



Here, we consider the simple test problem of a single particle interacting with a 
cloud of diffusing impurities, with the dynamics constrained in one dimension (a line). 
We restrict ourselves to the region of the parameter space in which the impurities have 
a vanishingly small probability to escape from the vicinity of the particle. Despite 
its apparent simplicity, such a system exhibits rich dynamics, but at the same time 
has features that make the problem analytically tractable. In the absence of external 
forces, we consider the statistics of monotonic excursions Ax of the particle, i.e. the 
distances the particle moves to a particular direction (here "left" or "right" along 
the one-dimensional line) without changing direction, see also figure HJ We find that 
these obey power law distributions P(Ax) ~ (Ax)~ TAx f c (Ax / Axq) with the exponent 
t~Ax = 4/3. The same is true for the monotonic changes of the external force AF when 
the particle is driven with a slow constant velocity. This paper is organized as follows: In 
the next section, we consider the interaction of a particle with a single mobile impurity, 
in the absence of external forces. Then we generalize this to the case with more impurity 
particles. In section [HI the effect of external drive is studied. Finally, section [4] finishes 
the paper with conclusions. 

2. A particle interacting with mobile impurities 

2.1. Single impurity 

As a starting point of our analysis, we consider the dynamics of a single particle 
interacting with one diffusing impurity particle. The equations of motion for the system 
are 



where x and x s are the positions of the particle and the impurity particle, respectively. 
f(z) is the interaction force between the particle and the impurity particle, \i defines 



types 0El[Tia. 



/id t x = f(x - Xg) 




(1) 
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the relative mobilities of the impurity and the particle and 77 is Gaussian white noise 
with standard deviation 5r] mimicking the effect of temperature. 

The dynamics of the particle can be analyzed by considering the stochastic process 
for the velocity d t x. By differentiating the equation of motion of the particle with respect 
to time and using the equation of motion for x s , one obtains 

li%x = d g f{z)[(l + t i)d t x-rj\. (2) 

Close to z = x — x s = 0, the force f(z) can be taken to be linear in z and thus the 
derivative of the force can be approximated by a constant, d z f(z) ~ — C, with C > 0. 
With A = C(l and £ = — (C//x)r/, equation (T5|) can then be rewritten in the form 

of an Ornstein-Uhlenbeck process for dtx, 

d 2 t x = -Xd t x + f . (3) 

For the process dtx, equation (jHJ) describes Brownian motion pushed toward the origin 
by a linear damping term. This problem has been considered e.g. in [11] , and the 
scaling exponents are known. In particular, the first return times T to origin of d t x 
have a probability distribution scaling as 

P(T) ~ T-^f c (p\ , (4) 



vT / 

with tt = 3/2 and the cut-off scale To ~ 1/A. Similarly, the average shape of an 
excursion, (d t x(t))T, scales for t,T — t <C 1/A as (d t x(t))T = T 7 " 1 f s ha P e(t/T), with 
7 = 3/2. Then, by using the scaling relation 7 = (tt — 1)/(tax — 1) [12], one obtains 
the distribution of the lengths Ax = J T d t xdt of monotonic excursions of the particle, 

P(Ax) ~ (Ax)-^f c , (5) 



VAxo 

with taz = 4/3. The cut-off scaling can be found as follows: The cut-off of the first 
return time distribution of equation (j3J) is given by T ~ 1/A [TT]. Here, both A and £ 
in equation ([3]) depend on C and /i. Due to the relation Axo ~ Tq with 7 = 3/2, Axo 
is expected to scale like 

We check this result in numerical simulations, where we take for simplicity f(z) 
to be of the form f(z) = — Az exp (— (l/2)(z//) 2 ), corresponding to d z f(z)\ z=0 = —A, 
i.e. C = A. We integrate the equations of motion (TIJ with the Euler algorithm. The 
strength of the thermal noise was chosen to be sufficiently weak such that the impurity 
cannot escape from the neighborhood of the particle. Figure [1] shows an example of the 
trajectories of the particle and the impurity. Figure [2] displays the distribution of Ax for 
various values of A. For a weak enough interaction strength A, the distributions display 
the expected scaling with taz = 4/3. The cut-off is found to scale as Aa;o ~ v4~ 1//2 , in 
agreement with our results above. 

This implies that the motion of the particle is reminiscent of truncated Levy flight, 
with a step length distribution given by equation ffih. However, the steps are not 
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Figure 1. An example of the trajectories of the particle and a single impurity, in 
the absence of external forces. Notice that the monotonic excursions of the particle 
correspond to motion of the particle during time intervals the impurity spends on a 
given side of the particle. Parameters of the simulation: \i = A = I = 1.0, 5r) = 0.1. 
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Figure 2. Main Figure: The probability distributions of monotonic excursions of the 
particle interacting with a single impurity, for various strengths A of the interaction 
force. For sufficiently low A values, the distribution scales with the exponent t& x — 
4/3, indicated by the solid line. Inset: A scaling plot of the distributions. The cut-off is 
observed to scale as Axo ~ y/l/A, in agreement with equation ((6|). Other parameters: 
1 = 1, (i=l and Si] — 0.1. 

instantaneous - their durations T exhibit power law scaling as well, equation (111) . By 
interpreting these step durations as waiting times between instantaneous steps, one 
would obtain for early times the scaling (x 2 )(t) ~ t 2 ( r T-i)/(TAa;-i) _ ^27 ^ However, 
here Ax and T are not independent (due to the relation (Ax) ~ T 7 ). As the early time 
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Figure 3. The mean square displacement of the particle interacting with a single 
impurity in the absence of external forces, for different strengths A of the particle- 
impurity interaction. The early times are characterized by superballistic motion with 
{x 2 )(t) ~ i 5 / 2 (solid line), while for longer times a cross-over to diffusive dynamics is 
observed (dashed line). Parameters of the simulation: fj, = I = 1.0, Sr) = 0.1. 



behavior is dominated by a single large step, one must consider instead the effect of a 
single step given its duration, 

(x 2 )(t)~ (\/\x(T)fP(T)dT = f T 2 ^ TT dT ~t 2 ~<- TT+ \ (7) 
Jo Jo 

corresponding to (x 2 )(t) ~ t 5//2 . Due to the truncated nature of the Levy flight, one 

expects a cross-over to diffusive behavior with (x 2 )(t) ~ t for long times. This is verified 

in figure El 



2.2. Several impurities 

The case of a fixed number iV > 1 of impurity particles is a straightforward 
generalization of that presented in the previous subsection. The equations of motion 
become 

fid t x = ^2f(x- x S)i ) 

i 

d t x S:i = - f(x - x Sti ) + 7]i. (8) 

From these equations, with the same procedure as above, one obtains an equation of 
the form of equation Q by setting A = C(N + n)/ fj, and £ = (C j fJ.)J2iVi- Thus, the 
same scaling, i.e. P(Ax) ~ (Arc) -4 / 3 is expected. One should notice, however, that 
the cut-off scale Axo is getting smaller with increasing N. With similar arguments as 
above, one finds that 
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Figure 4. Main Figure: The probability distributions of monotonic excursions of 
the particle interacting with N impurities. The distribution scales with the exponent 
t~Ax = 4/3, indicated by the solid line. Inset: A scaling plot of the distributions. The 
cut-off is observed to scale as Axo ~ 1/iV, in agreement with equation |9]). Parameters 
of the simulation were A = 0.01,Z = 1, fi = 1 and Si] = 0.1. 



This is verified by the numerical results presented in figure HI 
3. Constant velocity drive 

Next we proceed to study the effect of a weak external force F on the dynamics of 
the particle. In this context a constant velocity drive is perhaps the more interesting 
form of driving as a small constant force with d t F = does not change the dynamics 
from the non-driven case: equation (j2J) remains the same even if a constant force term 
is introduced in equation (pQ). In particular, we consider a particle driven by a force 
given by F = K(Vt — x), where V is the driving velocity and K is a spring constant 
characterizing the response of the driving mechanism. The equations of motion read 

fxd t x = J2f(x- x s>i ) + F 

i 

d t x S}i = - f(x - x s<i ) + r)i. (10) 

In systems like this driven with a constant velocity, the interesting quantity is the 
statistics of the external force fluctuations. To this end, we consider the stochastic 
process d t F. With a similar approach as above, one can write 



d 2 t F= -Kd 2 t x 



- + -(iV + M) 
fj, fx 



<9^ + — 5> 

^ i 



+ ^[V(N + ri-F], (11) 
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Figure 5. An example of the behaviour of the force as a function of time. The inset 
shows a magnification of a part of the signal in the steady state. Parameters of the 
simulation: fi = I = 1.0, A = 0.01, N = 10, K = 0.1, 6rj = 0.1 and V = 0.001. 



where the relation dtx = V — dtF / K has been used. In the steady state, the last term 
in equation (TTTT) has a zero mean, as one can write for the average steady state force 
F s = (fid t x — J2i f{% — x s,i)) = t^V + Nf s , where f s is the magnitude of the average 
retarding force acting on the particle due to a single impurity In the steady state the 
condition (d t x Sj i) = f s = V holds, and thus {V(N + ji) — F) =0. Assuming that the 
fluctuations <5[V(iV + /i) — F] = 5F are small compared to those of the white noise term 
in equation (TTTT) . i.e. 5F <^ VN5t], equation (flT|) can be approximately written as 

KC 

d t F + 5>, (12) 



d'fF 



- + -(iV + /i) 



which is again of the same form as equation ([3]). Thus, the monotonic changes of the 
external force AF = J T d t Fdt are expected to be distributed according to a power law 
P(AF) = (AF)- TAF f c (AF/AF ), with the exponent t af = 4/3 and the cut-off scale 
AF scaling as 

KCM (13) 
[K + C(N + fi)\ ^ 



Notice that the condition SF -C v NSr] implies that 

! ^VE « ! (14) 

[K + C(N + f,)] 3 ' 2 

For most of the relevant parameter values condition (|14p is fulfilled, only for KC ^> 1 
this is not the case. 

Next we check these predictions numerically for different values of the parameters 
satisfying the condition (fl4l) . Figure [5] shows an example of the behaviour of the force 
as a function of time. After an initial transient, the system reaches a steady state in 
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Figure 6. Distributions of the monotonic changes of the external force AF (upper 
panel) for various values of the parameters. The lower panel shows a data collapse 
of the distributions, with AFq computed from equation (fl3|) . All other cases have 
8r] = 0.1, while the data shown with black stars has 5i] = 1.0. Other parameters of the 
simulations: K = 0.1, A = 0.01 and V = 0.001. 

which the force fluctuates around a constant average value. In figure El we show the 
distributions of the monotonic changes of the external force AF in the steady state, 
for different values of the various parameters. Power law scaling of the distributions 
consistent with the exponent value taf = 4/3 is observed, with a cut-off of the 
distributions in agreement with equation (jTBl) . The different AF- values appear to be 
uncorrelated in time. 

4. Conclusions 

In this paper we have studied the dynamics of a single particle interacting with a 
cloud of diffusing impurities. In the absence of external forces, the problem can be 
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mapped to Brownian motion in a potential within the harmonic approximation of the 
attractive particle-impurity interaction. Thus, the monotonic excursions of the particle 
are distributed as a power law, P(Ax) ~ (Ax)~ TAx , with tax = 4/3. For a particle driven 
with a small constant velocity (such that the particle is dragging the impurity cloud 
without escaping from it), the external force fluctuations follow the same dynamics, 
which makes it possible to derive the probability distribution for the monotonic changes 
of the external force scaling as P(AF) ~ (AF)~ TAF , again with t^f = 4/3. 

While a typical experimentally relevant scenario would correspond to either a higher 
dimensional object such as a flexible line or a large number of interacting particles 
interacting with mobile impurities, the simple setup of the present study serves as a 
convenient starting point for such considerations illustrating the relevant phenomena 
in a transparent manner. One interesting observation is that already at the level of a 
single particle interacting with one or more mobile impurities the dynamics has scale 
free features, arising from the properties of simple random walks. Physical situations in 
which these kind of considerations could be relevant include dislocations interacting with 
solute atoms, a subject that has recently attracted considerable attention [13], [HI [15] . In 
many of these studies, more realistic interaction forces between dislocations and solute 
atoms have been used, but the focus has been on different quantities such as the average 
velocity of the dislocation. 

The most intriguing phenomena associated with the presence of diffusing impurities 
are the collective effects arising from the simultaneous interaction of large number of 
entities with each other and with mobile impurities. An example of such a system is 
provided by interacting dislocation ensembles interacting with diffusing solute atoms in 
solid solutions, giving rise to phenomena such as the Portevin-Le Chatelier effect. In the 
PLC effect, large numbers of dislocations synchronize their motion to form macroscopic 
deformation bands of various kinds. As this is widely believed to be due to the dynamic 
interaction of the dislocations with the diffusing solute atoms, a natural future line 
of research would be to study such effects by considering numerically the dynamics 
of a large number of interacting dislocations interacting with diffusing impurities. One 
further motivation for such studies could be the recent observation that purely stochastic 
effects can induce switching between collective motion states [16] . 
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